﻿ Harry Van Trees Kristine L Bell Steven Kay DETECȚIA SEMNALELOR DETERMINISTE CUNOSCUTE Vom aborda problema detecției unui semnal cunoscut Un exemplu tipic pentru un astfel de caz, este transmisia digitală coerentă, aplicată în telecomunicații Testul, pe care se bazează detecția, va rezulta sub o formă liniară, cunoscută și sub denumirea de filtru adaptat la forma semnalului, sau, simplu, filtru adaptat Termenul corespunzător în engleză este cel de matched filter Efectul unui astfel de filtru este cel de maximizare a SNR În figură (wikipedia) se sugerază efectul filtrului asupra unui semnal afectat, aditiv, de un zgomot alb, gaussian (AWGN) Filtrul adaptat la forma semnalului (the matched filter) Vom începe construcția detectorului de tip filtru adaptat, aplicând abordarea Neyman- Pearson (NP) Și abordarea bayesiană conduce la aceeași formă a detectorului, diferind doar pragurile de comparare și, evident, performanțele statistice Problema constă în a alege între două ipoteze H 0: [ ] [ ] ; n=0, 1, , N-1xn wn= H 1: [ ] [ ] [ ]; n=0, 1, , N-1xn sn wn=+ 2 wnσ [ ] 0, ∼N () zgomotul w[n] fiind alb, gaussian, cu secvența de autocorelație de forma 1 2 δδ=+== {}() 0[] [][ } [] 2 []wrk Ewnwn k k N kσ Reamintim că 1 , nn=≥ 0 1 , 0 ⎧⎧ [ ] [ ]nunδ== ⎨⎨ 0 , 0 0 , 0nn≠ Lγx () () pxH 0; Datele x fiind gaussiene, e suficient să cunoaştem mediile şi dispersiile lor, în cele două ipoteze H {}{} 0: [ ] [ ] =0Exn Ewn= = + = H {}{ }{}{ } 1:[][][][] [][]E xn E sn wn E sn E wn sn=+ 2 = H {} { } 0:[] []=Disp x n Disp w nσ 2 =+= = H {} { } { } 1:[] [][] []Disp x n Disp s n w n Disp w nσ Repartiţiile datelor, x, în cele două ipoteze, rezultă că sunt 1N− ⎫ 211⎧ − xH ()() ⎨⎬ ∑ exp [ ] [ ]pxnsn=−122; πσ () = 022nσ ⎩⎭ 1N− 211⎧⎫ xH ()() ⎨⎬ ∑ exp [ ]pxn=−022; πσ () ⎩⎭ = 022nσ Dacă substituim acest ultim rezultat în raportul de plauzibilitate, condiţia de selectare a ipotezei H1devine 1H 11NN−− ⎡⎤ 221⎧⎫ >x =− − − ()()() ⎨⎬ ∑∑ 2exp [ ] [ ] [ ]Lxnsnxnγ ⎢⎥ == 002nnσ ⎣⎦ ⎩⎭ sau, după logaritmare 1H 11NN−− 21⎡⎤ >xx ==− − () ()() ∑∑ 2ln [ ] 2 [ ] [ ] lnl L sn xnsnγ ⎢⎥ == 002nnσ ⎣⎦ Inegalitatea se pune sub forma 1H 11NN−− 211 [][][] lnxnsnsnγ −> () ∑ 22∑ σ 002nnσ == Deoarece s[n] este un semnal cunoscut, ce nu depinde de datele x, termenul marcat poate fi încorporat în prag aşa că 1H 11NN−− 221 ′ [][] l[]nsnTxnsnσγγ=+=>x ()() ∑∑ 002nn== sau 1H 1N− ′ = [ ] [ ] (12)Txnsnγ>x () ∑ 0n= 2 Detectorul NP obţinut constă dintr-o statistică a testului, T(x), ce depinde de realizările datelor, x, şi dintr-un prag al testului, γ ’ Valoarea pragului trebuie aleasă astfel încât să fie satisfăcută constrângerea impusă de valoarea acceptată pentru probabilitatea alarmei false Vom considera un exemplu simplu, în care semnalul cunoscut este o componentă continuă, s[n]=A>0 Substituind în test obţinem 1N− [] , []snA T Axn=⇒ =x () ∑ 0n= Se substituie statistica în testul (12)şi împărţind cu NA>0, ceeace nu schimbă sensul inegalităţii, obţinem 1H 1N− 11 1 ′′′′ []TT xnγ γ>xx == = ()() ∑ 0nNA N NA= Observaţie ′′ : 0xγ >>H A>0, decidem 1Pentru dacă ′′ : 0xγ x () ∑ 0n= Datele x fiind gaussiene, T(x) are şi ea o repartiţie gaussiană Este deci suficient sa-i determinam mediile şi dispersiile, în ambele ipoteze Pentru medii avem relaţiile TTT ==sx sw s wHH {}{} {}{} ; 00; 0ET E E E== TTTTT +=+==sx s s w ss s w ssHH ε {}(){} {}{} ET E E E== ; 11; Putem determina acum şi dispersiile 2 TTT = sx sw swHH {} {}{}() Disp T Disp Disp E== ; 00; {} 2TT T T ε EE== = swws s ww sσ {}{} TT + sx s s wHH {}() {}{} Disp T Disp Disp== ; 11; 22TT ε Disp== = sw ssσσ {} În cele două ipoteze, statistica T(x) pe baza căreia se iau deciziile, are, în cele două ipoteze, distribuţii normale, adică 2 ⎧ NHε () σ , în ipoteza 00, ⎪ ∼ T⎨ 2 εε NH () ⎪ σ , în ipoteza 1, ⎩ Suntem în cazul în care cele două ipoteze se disting prin mediile repartițpplor Ca atare ⎛⎞ ′ γ ′ H => = {} FAPPT Qγ 0; 2⎜⎟ ε σ ⎝⎠ și ⎛⎞ ′ − γε ′ => =H {} DPPT Qγ 1; 2⎜⎟ σ ε ⎝⎠ Din prima ecuație se determină pragul testului ca fiind 21− ′ = ε () FAQPγσ Se substituie expresia pragului în expresia probabilității de detecție și obținem ⎛⎞ ε 1− () ⎜⎟ 2DFAPQQP=− σ ⎝⎠ Dacă semalul de detectat este o componentă continuă, s[n]=A, energia se calculează cu 2 NA=ε 6 În figură se prezintă curbele de performanță ale detectorului NP, cu filtru adaptat la forma semnalului Utilizarea filtrului adaptat ne conduce la conceptul de “câștig de prelucrare” El poate fi privit ca descriind avantajul obținut dacă deciziile se iau pe baza unei statistici, față de cazul în care deciziile se iau direct, pe baza datelor primare Câștigul de prelucrare se cuantifică prin creșterea SNR, obținută prin prelucrare La intrarea filtrului adaptat, SNR este 2 A = 2inη σ La ieșirea filtrului adaptat, ca urmare a prelucrării a N eșantioane, SNR devine 2 NA == max2outηη σ Măsurat în dB, câștigul de prelucrare (Processing Gain) este η 10log 10loginPG N== outη Prelucrând nu unul, ci 4 eșantioane, obținem un câștig de prelucrare de 4dB Reamintim că, în prezentarea filtrului adaptat, am presupus că zgomotul este alb, gaussian Se poate pune problema conceperii unui filtru adaptat pentru zgomot gaussian colorat, caracterizat de o matrice de covarianță C 7 Filtrul adaptat generalizat (rezumat) Dacă eșantioanele de zgomot sunt corelate, atunci w0C∼N, () Elementele matricei de covarianță au expresiile [][] [ ]Ewmwn r m n==−C {} [] wmn Drept urmare a formei funcției de corelație a zgomotului, matricea de covarianță are o formă Toeplitz, simetrică Vom presupune că urmărim să detectăm semnalul cunoscut, s, afectat, aditiv, de zgomotul colorat, w Datele recepționate au, în cele două ipoteze, repartițiile 0CNH () ⎧ , în ipoteza 0, x∼ ⎨ H sCN () ⎩ , în ipoteza 1, Densitățile de repartiție se determină cu relațiile ⎧⎫ 111T− xxsCxsH −− ()() () ⎨⎬ expNp=−12; 22π ⎩⎭ C () ⎧⎫ 111T− xxCxCH () ⎨⎬ 02;Np=− exp ; unde este determinantul 22π ⎩⎭ C () Pentru a lua o decizie, vom construi criteriul raportului de plauzibilitate logaritmică 1H () pxH 1; >γ ln lnlnlL==xx () () () pxH 0; Prin calcul direct se arată că raportul de plauzibilitate poate fi pus în forma 1H 111TT−− lnlγ=−>sC sxxCs () 2 Termenul marcat nu depinde de date, așa că poate fi inclus în prag 1H 111TT−− Tγγ=+=>xxCs sCsln (14) () ′ 2 Detectorul definit prin relația (14), se numește corelator de replică generalizat,sau filtru adaptat (la forma semnalului), generalizat Dacă transformăm mai întâi semnalul cunoscut, s, prin înmulțire cu inversa matricei de covarianță 1− ′ =sCs atunci corelatorul de replică generalizat poate fi considerat ca fiind un corelator de replică obișnuit, dar care are ca referință semnalul “distorsionat”, s’ 1TT− T′==xxCsxs () Pentru exemplificare vom considera că eșantioanele de zgomot deși nu sunt corelate, au dispersii diferite între ele Matricea de covarianțăși inversa ei sunt 8 2 ⎡⎤ 000σ 2⎢⎥ −− −− 22 2 2 2 21100σ ⎢⎥ == = CC…… σ σ σ σ {}{ } −− 01 1 0 1 1,,, ; , ,,NNdiag diagσσ ⎢⎥ 2⎢⎥ − 100Nσ ⎣⎦ Substituim în forma (14) a detectorului și, efectuând calculele obținem 2− ⎡⎤ ⎡⎤ 0 00sσ ⎢⎥ ⎢⎥ 2− − 11 00Tsσ ⎢⎥ [ 1]TxxxN⎢⎥ == − xxCs… () [] ⎢⎥ ⎢⎥ ⎢⎥ ⎢⎥ 2− − ⎣⎦ 1 00NsNσ− ⎣⎦ 1H 1N− [][]xnsn => ′ ∑ 2γ 0nnσ= *** Reamintim că se poate aplica tehnica de “albire” a datelor Avem 1T− C=DD Pentru exemplul anterior avem T DD 211−−− ⎡⎤⎡⎤⎡⎤ σ 00000 00 00σσ ⎢⎥⎢⎥⎢⎥ 211−− − σ − 1111000000σσ ⎢⎥⎢⎥⎢⎥ C == ⎢⎥⎢⎥⎢⎥ ⎢⎥⎢⎥⎢⎥ 211−− − σ 11100 00 00NNNσσ −−− ⎣⎦⎣⎦⎣⎦ Statistica testului se poate pune și sub forma 1H 1TTTT− ′′ Tγ>x x C s x D Ds Dx Ds x s == = = () ( )( ) ′ xDx sDs ′′ si == Reamintim că matricea D se numește matrice de albire Datele se decorelează prin înmultirea cu matricea de albire, dar și replica memorată a semnalului căutat trebuie înmulțită cu D În rest avem de-a face cu un corelator de replică obișnuit, sau cu un filtru FIR adaptat obișnuit 9 Performanțele statistice ale filtrului adaptat generalizat Se arată că statistica (14) 1H 1T− Tγ=>xxCs () ′ are, în cele două ipoteze, repartițiile 1T− sC sN H () ⎧ , în ipoteza 00, ⎪ T⎨∼ 11TT−− sC s sC sN H () ⎪ , în ipoteza 1, ⎩ Cele două repartiții diferă prin medie, așa că vom calcula probabilitatea de detecție servindu-ne de relația stabilită anterior 12− ( 15)PQQP d=− () ) ( DFA în care coeficientul de deflexie are epresia 2 − () 210μμ 2d= σ Substituind găsim 21T− d=sC s și 11T−− sC s () ) ( DFAPQQP=− În cazul în care zgomotul este alb, probabilitatea de detecție depinde de raportul dintre energia semnalului util și puterea zgomotului, dar nu și de forma semnalului În cazul zgomotului colorat, coeficientul de deflexie depinde de forma semnalului Pentru a obține o probabilitate de detecție mare, am putea maximiza d, fără a crește energia semnalului cu care se face transmisia, doar alegând o formă adecvată Este necesar să maximizăm expresia 21T− d=sC s cu constrângerea T =ssε Se costruiește lagrangeanul 1TT− Fλ=+−sC s ssε () Reamintim regulile de derivare TT bx xAx T∂∂ ==; 2 , dacă: =bAxAA ∂∂xx Cu acestea 1TT− ∂∂ ∂sC s ss 1F− 220λλ =−=−=Cs s ∂∂ ∂ss s sau, în final 1− λ=Cs s 10 ceeace înseamnă că s este un vector propriu al inversei matricei de covarianță iar λ este o valoare proprie atașată acelui vector propriu Din relația de constrângere rezultă 1TT− λλ==sC s ssε Membrul stâng trebuie maximizat; energia fiind constantă, maximul se obține pentru valoarea proprie, λ, maximă Ne putem referi la matricea de covarianță, nu la inversa ei În acel caz se arată că va trebui să alegem valoarea proprie, minimă, μ=1/λ Atragem atenția că vectorii proprii se determină cu norma unitară, motiv pentru care trebuie ponderată cu energia admisă pentru semnalul de transmisie Pentru exemplificare vom considera că eșantioanele de zgomot sunt corelate, având matricea de covarianță 1ρ ⎡⎤ = C ⎢⎥ 1ρ ⎣⎦ limitată la dimensiunile 2x2, deoarece la recepție se lucrează doar pe două eșantioane de date xx=x [] Se determină valorile proprii ale matricei de covarianță cu ecuația caracteristică 21μρ 2− ρμρμρ =− − = =+ =− () 1210; 1 1μ 1ρμ− Se determină cei doi vectori proprii de normă unitară cu relațiile − ⎡⎤⎡⎤ ⎡⎤ 12210iivμρ == 12; 1; 1, 2iivv i=+ ⎢⎥⎢⎥ ⎢⎥ − ⎣⎦ 210iivρμ ⎣⎦⎣⎦ Vectorii proprii și valorile proprii aferente sunt ⎡ ⎤ 12 12⎡⎤ ρμρ==+ = =− ⎢⎥ ⎢ ⎥ 112 2; 1 si ; 1vvμ 12 12− ⎢⎥ ⎢ ⎥ ⎣⎦ ⎣ ⎦ Vom presupune că avem ρ>0, (dar se poate considerași cazul opus) Valoarea proprie minimă și vectorul propriu corespunzător, precum și cele două eșantioane cu care se recomandă să efectuăm transmisiile sunt, în acest caz 12⎡⎤ 01ρμ μρε >⇒ = =− ⇒=s ⎢⎥ Dacămin 2 12− ⎢⎥ ⎣⎦ Cele două esantioane ale semnalului su 2; 2;ntss== − εε Statistica testului care maximizează probabilitatea de detecție devine ⎡⎤ 11211T− = = −xxCsxεε 2 Txx xx== ()() [] ⎢⎥ − − min112μρ ⎢⎥ ⎣⎦ 11 Decizia se va lua comparând o diferență ponderată a celor două eșantioane de date cu un prag Din cauza corelării pozitive a eșantioanelor de zgomot, ele vor avea, statistic, același semn și valori apropiate În consecință, efectul zgomotului asupra statisticii T se reduce, nu însă și efectul semnaluilui util, după cum rezultă din diferența 2 2+ xx sw sw ww−= + − + = −ε ()()() μ deflexia Cu cele stabilite, se poate calcula coeficientul de deflexie Deoarece λ=1/ devine 211T− sCsε ;d== − 1ρ s maxλ minμ valoare ce trebuie utilizată în relația (15) pentru calculul probabilității de detecție, la o valoare impusă a probabilității alarmei false Este de remarcat că 2 pentru1, dρ→→∞ și, în consecință 1P → D Explicația este simplă Cu cât corelarea eșantioanelor succesive de zgomot este mai apropiată de 1, cu atât mai mult scade, statistic, diferența lor, ceece duce la creșterea performanțelor de detecție 12